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The effect of an external magnetic field on the finite temperature transition of QCD is studied. We 
measure thermodynamic observables including the quark condensates and susceptibilities and 
the strange quark number susceptibility. We generate configurations at various values of the 
quantized magnetic flux with Nf = 2 + 1 flavors of stout smeared staggered quarks at physical 
quark masses. We perform the renormalization of our observables and approach the continuum 
limit with N t =6,8 and 10 lattices. We also check finite volume effects using various lattice 
volumes. Our main result is that the transition temperature significantly decreases with growing 
magnetic field, and that the transition remains an analytic crossover up to our largest external field 
JeH w 1 GeV. 
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1. Introduction 

The response of the QCD vacuum to external (electro)magnetic fields is relevant for several phys- 
ical situations, since very strong magnetic fields are thought to be present during the electroweak 
phase transition in the early universe |jl|], in the interior of dense neutron stars called magnetars 
and in non-central heavy ion collisions [^]. In the latter case, although the external field has a very 
short lifetime (of the order of 1 fm/c), the magnetic 'impulse' coincides with the generation of the 
quark-gluon plasma and thus may have a significant effect on the properties of the transition. 

The structure of the phase diagram of QCD in the magnetic field-temperature (B-T) plane has 
been studied extensively in the past years. Calculations have been carried out within various low 
energy effective models of QCD. Most of these models predict an increasing transition temperature 
T c and a strengthening of the transition with growing B, see, e.g. [|L BL [fl. However, the opposite 
effect of a decreasing deconfinement transition temperature was predicted using chiral perturbation 
theory for two quark flavors Moreover, in most of these model calculations chiral symmetry 



breaking is enhanced with growing B through an increase in the chiral condensate \^ |lfj, |11|]. On 
the other hand, it was also conjectured that the running of the strong coupling in the presence of 
magnetic fields may modify this magnetic catalysis, and even turn the effect around to make the 
dynamical mass decrease with B in some regions Ji2|]. 



In a recent lattice simulation with Nf = 2 flavors of staggered quarks [ ]13| ] the chiral conden- 
sate was observed to grow with the external field for any temperature T in the transition region. 
The size of this effect was however found to be different for different values of T, resulting in 
an increase in both the pseudocritical temperature T c and the strength of the transition, being in 
qualitative agreement with most of the model predictions. In this project our aim is to perform a 
similar lattice study, but with improved gauge and smeared fermionic actions and with Nf = 2 + \ 
flavors of quarks, at the physical pion mass, and extrapolate the results to the continuum limit. We 
explore a wide temperature region around the zero-field pseudocritical temperature T C (B = 0), for 
various values of the magnetic field, ranging up to \fe~B ~ 1 GeV, i.e. covering the region that 
is phenomenologically interesting for noncentral heavy ion collisions and for the evolution of the 
early universe. Our results are also published in a separate paper 

2. Magnetic field on the lattice and observables 

We consider a constant external magnetic field B = (0,0, B) in the z direction, which can be re- 
alized by a vector potential A y = Bx in the continuum. On the lattice such a vector potential can 
be represented by complex phases u v (n) € U(l) that multiply the U v (n) G SU(3) links of the lat- 
tice d|, n, jni, u, Hi, 



u y (n)=e ia « Bn *, u x (N x -l,n y ,n z ,n t )=e- ia « BN < n >, 

\£. i ) 

u x (n) = 1 for n x ^ N x — 1, u v {n) = 1 for v {*,;y}, 

where the sites are labeled by integers n = (n x ,n y ,n z ,n t ), with n v = . . .N v — I, and A^ v is the num- 
ber of lattice points in the direction v. In this formulation we have periodic boundary conditions 
in all spatial directions and the magnetic flux going through any plaquette in the x — y plane is 
constant. 
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It is well known that in a finite box with periodic boundary conditions the magnetic flux is 
quantized in terms of the area of the plane orthogonal to the external field [19, 20]. On the lattice 
this leads to the quantization condition, 

aV? = _±, N b €Z, (2.2) 

where q is the smallest charge in the system - in our case the down quark charge q = q#. We 
note that the lattice magnetic field is periodic in Nb with a period of N x N y . In order to have an 
unambiguous implementation we constrain the flux such that < Nb < N x N y /4. At larger values 
of Nb the periodicity is expected to introduce saturation effects, like it was observed in [21, 22]. 

We derive our observables from the staggered partition function with three flavors u,d and 
s. The quark flavors are treated separately since their charges/masses are different: q u = —Iqj = 
—2q s , and we assume m u = mdi^ m s . The dependence on the external field enters the partition 
function only through the determinants, in the form qfB. Since we are concerned with a constant 
external field we do not include any dynamics for the U(l) field introduced above. However our 
approach is 'dynamical' (i.e. not quenched) in the sense that the magnetic field is taken into account 
both in the configuration production and in the measurements. 

Our observables include the chiral condensates and susceptibilities for the light flavors / = 
ii. d, and the strange quark number susceptibility, 

where we defined the spatial volume of the system as V = (N s a) 3 with N s = N x = N y = N z . The 
condensate will be denoted in the following by the first letter of the flavor name, e.g. uu. 

In order to approach the continuum limit, the renormalization of these observables has to be 
carried out. We remark that the divergences to be canceled are independent of both T and B [23, 13]. 
Therefore we can eliminate the additive divergences by subtracting the T = 0, B = contribution. 
Then we multiply by an appropriate power of the bare quark mass to cancel the multiplicative 
divergences 



Y ¥ f(B,T)=m f [xfr ¥f (B,T)-rf f¥f (0,0)] -1 X f(B,T) = m) [Xf(B,T) -#(0,0)] -\. (2.4) 

This procedure leads to a renormalized condensate that, for B = 0, is zero at T = and approaches 
a negative value as T is increased. Considering the strange quark number susceptibility, c s 2 needs 
no renormalization (neither at B = nor at B ^ 0) since it is connected to a conserved current. 

In this study we use the free-level improved Symanzik gauge action and stout smeared stag- 
gered fermions; details about the action can be found in [p5|]. We generate lattice configurations 
both at T = and T > with an exact RHMC algorithm, for various values of the gauge coupling 
and the magnetic flux. For our finite temperature runs we have lattice configurations with N t =6,8 
and 10. Finite volume effects are studied on the N t = 6 ensemble using sets of N s = 16, 24 and 32 
lattices. The masses of the up, down and strange quarks are set to their physical values along the 
line of constant physics (LCP) by fixing the ratios //f/m^ and /k^/mjc to their experimental values. 
The lattice spacing is determined by fy. Details of the determination of the LCP and the lattice 
scale can be found in, e.g. [26]. Since the finite external field does not affect the lattice scale [|T4|], 
we use the lattice spacing measurements at T = and B = to set the scale also at T ^ and B^0. 
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3. Chiral condensate at finite T and B 



0.1 
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In figure [j] we plot the renormalized chiral conden- 
sate as a function of the temperature for three values 
of B. The transition temperature apparently decreases 
with increasing B, thereby contradicting a vast number 
of model calculations, see the summary given in the in- 
troduction. Furthermore this observation also disagrees 
with the lattice result of [JT3Q. First of all, to check our 
simulation code we reproduced the results of Jl3| ] with 
an identical setting (i.e. same lattice sizes and quark 
masses, no smearing and no improvement in the ac- 
tion) at a couple of temperatures. Since we find a per- 
fect agreement, we conclude that we are left with three 
possible reasons for the discrepancy. First, the lattice 
spacing of [|13|] is larger, a « 0.3 fm, and also an unim- 
proved action is used, so lattice discretization errors may be significant. Second, the present study 
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Figure 1: The renormalized chiral conden- 
sate as a function of T for three different 
values of B for the N t =6 lattices. The tran- 
sition temperature shifts towards the left. 



uses Nf ■ = 2 + 1 flavors as opposed to the Nf = 2 of Q13Q, and the pseudocritical temperature is 



known to depend on the number of flavors [|27|], which may also introduce systematic differences 



in the dependence on the external field. Third, the quark masses of [ ]13[ ] are larger than in the 
present study, which can also cause drastic changes in thermodynamics - for example the nature of 
the transition at B = depends very strongly (and non-monotonically) on the quark masses. 



On closer inspection, the differences between our results and those of []13j] can actually be 
traced back to the behavior of the chiral condensate as a function of B and T. While the authors 



of [ ]13[ ] observed that at any temperature the condensate increases with B, we find that this depen- 
dence is more complex and in the transition region the condensate decreases with growing B. As 
already mentioned in the introduction, the possibility of such a decrease in the condensate with B 
was also raised in low energy model calculations [|12|]. We find that for three strange quarks this 
complex dependence turns into a monotically increasing behavior for all temperatures Jl4|], which 
indicates that the response of the chiral condensate to the external field is very strongly influenced 
by the quark masses. We summarize our findings as a) the dependence of the condensate on the 
external field is non-monotonic and varies strongly with temperature, and b) as a result the pseu- 
docritical temperature shifts to lower values at large B as compared to the B = case. The latter 
observation is supported by a similar T C (B) dependence deduced from the chiral susceptibility, the 



Polyakov loop and the strange susceptibility [14]. 



4. Transition strength and transition temperature 

First we study the strength of the transition as a function of the external field. At B = the transition 



is known to be a broad crossover [28], where the approximate order parameters like the chiral 
condensate and the Polyakov loop change smoothly with the temperature, and no finite volume 
scaling is visible in the observables. 
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Figure 2: The chiral susceptibility as a function of T measured on our N t = 6 lattices for different spatial 
volumes (left panel). Relative changes in the T-dependence of %u as measured on Nt = 6 lattices (right 
panel). 

To properly determine the nature of the transition we search for finite volume scaling in our 
observables. To this end we perform simulations at our largest magnetic field on the N t = 6 lattices 
with /V v = 16,24 and 32. The largest lattice in the transition region corresponds to a box of linear 
size ~ 7 fm. Here we keep eB/T 2 fixed (and not B itself) as we are only interested in differences 
between the various volumes. In the left panel of figure || the results for the chiral susceptibility are 
shown as a function of the temperature for eB/T 2 82. The figure shows that our = 16 results 
agree within statistical errors with the N s = 24 and N s = 32 data, indicating that finite size errors 
are small, compared to statistical errors. This independence on the volume also implies that the 
transition at this high magnetic field is still an analytic crossover. 

To further study how the strength of the transition changes we investigate the width of the 
renormalized chiral susceptibility. In the right panel of figure |2] we plot the susceptibility divided by 
its maximum value as a function of T — T C (B) for three different values of the magnetic field for the 
N t = 6 lattices. We find that the width of the peak is only mildly affected by the magnetic field. In 
particular, the width of the peak at half maximum decreases from ~ 30(3) MeV to ~ 25(3) MeV as 
the external field is increased from zero to eB = 1.05 GeV 2 . We find a very similar behavior on the 
N t = 8 and 10 lattices. From this analysis our final conclusions are that the width of the transition 
decreases only mildly with increasing magnetic field, and as the finite size scaling analysis has 
shown, the transition remains an analytic crossover at least up to \fe~B ~ 1 GeV. 

Next we study the observables as functions of the temperature to determine the pseudocritical 
temperature. We search for the inflection point of the renormalized chiral condensate uu r + dd r 
and the strange quark number susceptibility c|. To carry out the continuum extrapolation, we fit 
the results for T C (B) for all three lattice spacings (N t = 6, 8 and 10) together with an N t -dependent 
polynomial function of order four of the form T c (B,N t ) = L,- = o( a i +biN t ~ 2 )B'. This ensures the 
scaling of the final results with Nf ~ a 2 . We obtain # 2 /dof. w 0.5... 1.2 indicating good fit 
qualities. The results are plotted for the condensate and the strange susceptibility in figure [3| 

As is clearly visible in the figure, both observables show that the pseudocritical temperature 
decreases with growing external field B. We mention that the chiral susceptibility gives a similar 
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phase diagram, and that preliminary results about the Polyakov loop also indicate a decrease in 
T e (B) @. 



5. Summary 

In this paper we studied the finite temperature transition of QCD in the presence of external (elec- 
tromagnetic fields via lattice simulations at physical quark masses. The extrapolation to the con- 
tinuum limit is carried out, and finite size effects are under control. The results are relevant for the 
description of both the evolution of the early universe and of noncentral heavy ion collisions. 

We obtained the phase diagram of QCD in the B — T plane in the phenomenologically in- 
teresting region of < eB < 1 GeV 2 . Performing a finite volume scaling study we found that 
the transition remains an analytic crossover up to our largest magnetic fields, with the transition 
width decreasing only mildly. This rules out the existence of a critical endpoint in the B — T phase 
diagram below eB = 1 GeV 2 . 

Moreover, our results indicate that the transition temperature significantly decreases with in- 
creasing B. This result contradicts several model calculations present in the literature which predict 
an increase in T c as B grows (see the summary in section [I]). By comparing our magnetic fields to 
the maximal fields that may be produced in noncentral heavy ion collisions we conclude that the de- 
crease in T c is negligible for RHIC and may be up to 5 — 10 MeV for the LHC. Moreover, the effect 
grows with the magnetic field, exceeding 20% for <f 2 at y/eB = 1 GeV. This may have a significant 
impact on the description of the QCD transition during the evolution of the early universe. 

Acknowledgments 

This work has been supported by DFG grants SFB-TR 55, FO 502/1-2 and BR 2872/4-2, the EU 
grants (FP7/2007-2013)/ERC no 208740 and PITN-GA-2009-238353 (ITN STRONGnet). Com- 
putations were carried out on the GPU cluster at the Eotvos University in Budapest and on the 
Bluegene/P at FZ Jiilich. We thank Ferenc Niedermayer for useful discussions, interesting ideas 
and for careful reading of the manuscript. G.E. would like to thank Massimo D'Elia, Swagato 
Mukherjee, Daniel Nogradi, Tamas Kovacs and Igor Shovkovy for useful discussions. 



6 



The finite temperature QCD transition in external magnetic fields 



G. Endrodi 



References 



T. Vachaspati Phys. Lett. B265 (1991) 258-261. 

R. C. Duncan and C. Thompson Astrophys. J. 392 (1992) L9. 

V. Skokov, A. Y. Illarionov, and V. Toneev Int. J. Mod. Phys. A24 (2009) 5925-5932, 
[ |arXiv: 0907 . 13 9 6p . 

A. J. Mizher, M. N. Chernodub, and E. S. Fraga Phys. Rev. D82 (2010) 105016, 
[larXiv: 1004 .27121]. 



E. S. Fraga and A. J. Mizher Nucl. Phys. A820 (2009) 103c-106c, [arXiv: 0810 . 3693], 



R. Gatto and M. Ruggieri Phys. Rev. D83 (201 1) 034016, [arXiv: 1012 . 1291], 



K. Kashiwa/%5. Rev. D83 (201 1) 1 17901, [arXiv: 1104 5167], 



N. O. Agasian and S. M. Fedorov Phys. Lett. B663 (2008) 445^149, [ |arXiv: 0803 . 315 6p . 

V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy Nucl. Phys. B462 (1996) 249-290, 
Ohep-ph/9 5 9 3 2 00. 



S.-i. Nam and C.-W. Kao Phys. Rev. D83 (201 1) 096009, [ arXiv : 1103 . 6057 |. 



J. K. Boomsma and D. Boer Phys. Rev. D81 (2010) 074005, [ arXiv : 911 . 2 1 64 |. 



V. A. Miransky and I. A. Shovkovy Phys. Rev. D66 (2002) 045006, [hep-ph/0205348]. 



M. D'Elia, S. Mukherjee, and F. Sanfilippo Phys. Rev. D82 (2010) 051501, [ |arXiv: 1005 . 5 3 65p . 
G. S. Bali, F. Bruckmann, G. Endrodi, Z. Fodor, S. D. Katz, S. Krieg, A. Schafer and K. K. Szabo 



arXiv: 1111.4956 



G. Martinelli, G. Parisi, R. Petronzio, and F. Rapuano Phys. Lett. B116 (1982) 434. 

C. W. Bernard, T. Draper, K. Olynyk, and M. Rushton Phys. Rev. Lett. 49 (1982) 1076. 

L. Zhou, F. X. Lee, W. W ilcox, and J. C. Christensen Nucl. Phys. Proc. Suppl. 119 (2003) 272-274, 
[ |hep-lat/0209128[ |. 

D. S. Roberts, P. O. Bowman, W. Kamleh, and D. B. Leinweber Phys. Rev. D83 (201 1) 094504, 
[ |arXiv: 1011 . 1975[ |. 

G. 'tHooft Nucl. Phys. B153 (1979) 141. 



M. H. Al-Hashimi and U. J. Wiese Annals Phys. 324 (2009) 343-360, [arXiv: 0807 .0630], 



M. D'Elia and F. Negro Phys. Rev. D83 (201 1) 1 14028, [arXiv: 1103 . 2080], 



F. Bruckmann and G. Endrodi Phys.Rev. D84 (201 1) 074506, [ |arXiv: 1104 . 5664[ |. 
A. Salam and J. Strathdee Nuclear Physics B 90 (1975) 203-220. 

G. Endrodi, Z. Fodor, S. D. Katz, and K. K. Szabo JHEP 04 (201 1) 001, 



arXiv: 1102 . 135< 



Y. Aoki, Z. Fodor, S. D. Katz, and K. K. Szabo JHEP 01 (2006) 089, [hep-lat/0510084]. 



S. Borsanyi et. al. JHEP 11 (2010) 077, [arXiv: 1007 . 2580], 



F. Karsch, E. Laermann, and A. Peikert Nucl. Phys. B605 (2001) 579-599, [ |hep-lat/0012023| ]. 

Y. Aoki, G. Endrodi, Z. Fodor, S. D. Katz, and K. K. Szabo Nature 443 (2006) 675-678, 
[|hep-lat/0 611014|]. 



7 



